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ABSTRACT 

We report on the detection of a molecular band centered at ~98 /im 
(~102 cm~^), observed with the Infrared Space Observatory^ in the young Plane- 
tary Nebula NGC 7027. The band structure and intensity can not be reproduced 
by atomic fine structure lines, recombination lines or by the rotational emission 
of abundant molecules. We discuss the possible contribution of the low-energy 
bending modes of pure carbon chains to the unidentified far-IR bands observed in 
C-rich evolved objects. In particular, we speculate that the band emission could 
arise from the z/g and z/7 bending modes of Cg and C5, for which wavenumbers of 
90±50 and 107±5 cm~^ have been estimated from photoelectron spectroscopy. 

Subject headings: infrared: stars — circumstellar matter — stars: individual 
(NGC 7027) — ISM: molecules — fine: identification 



1. Introduction 



Carbon chains molecules are of special interest for Astrophysics since Douglas (1977) 
proposed them as possible carriers of the diffuse interstellar bands (DIBs). In the past 
years, many polar carbon chains have been detected in the interstellar and circumstellar 
media (ISM and CSM) through their pure rotational spectrum at radio wavelengths. Major 
examples are the cyanopolyynes, HC2n+iN with n=l-5 (Turner 1971, Avery et al. 1976, 
Kroto et al. 1978, Broten et al. 1978, Bell et al. 1997), and the hydrogenated carbon chain 



^Based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe- 
ciahy the PI countries: France, Germany, the Netherlands and the United Kingdom) and with participation 
of ISAS and NASA. 
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radicals such as C5H, CeH, C7H and CsH (Cernicharo et al. 1986, Guelin et al. 1987, 1997, 
Cernicharo & Guelin 1996). 

On the other hand, larger carbon molecular complexes such as the polyatomic aromatic 
hydrocarbons (PAHs) are believed to dominate the ubiquitous mid-IR emission seen in the 
unidentified infrared bands (UIBs). However, only the single aromatic species, benzene, has 
been so far identified through its IR active modes (Cernicharo et al. 2001). 

All this chemistry richness suggests that the growth mechanisms producing complex 
carbon molecules such as the PAHs or the fuUerenes are highly efficient. Nevertheless, the 
set of 'building blocks" and the possible chemical reactions still have to be identified and 
clarified. Among the possible species, the pure carbon chains, (C^; Van Orden & Saykally 
1998), could be the 'skeletons" from which larger organic molecules can be formed (see 
Cernicharo, Goicoechea & Caux 2000, Cernicharo 2004). Due to the lack of permanent 
electric dipole, these species do not have rotational spectrum to be observed from radio 
telescopes. The only way to detect them in the dense ISM and GSM is through their 
asymmetrical stretching modes around ~5 /im (2000 cm~^) and/or through their low-energy 
bending modes around ~100 //m (~100 cm~^). The first technique allowed the detection of 
C3 (Hinkle et al. 1988) and C5 (Bernath et al. 1989) in the circumstellar envelope of the 
bright IR evolved star IRC+10216. However, few sources (specially in the ISM) have enough 
flux at ~5 /im to allow systematic studies of the C„ stretching modes, and the less known 
far-IR bending modes are the only way to detect these species. 

Before the launch of ISO (Kessler et al. 1996), we proposed to observe the bending 

modes of several polyatomic molecules in the far-IR. As a result, the z/2 bending mode of 
C3 has been observed in Sgr B2 and IRC+10216 (Cernicharo et al. 2000). In addition, 
we have detected several Unidentified far-IR bands (UflBs) in the ISO spectrum of many 
C-rich sources that may also be related with the bending modes of carbon chains species. 
In particular we have tentatively assigned a UfIB at ~57.5 fj,m (~174 cm~^), observed in 
Sgr B2 and in C-rich evolved stars including the planetary nebula (PNe) NGG 7027, to the 
i/5 bending mode of G4 (Cernicharo, Goicoechea & Benilan 2002). 

NGC 7027 is one of the most studied evolved objects in the Galaxy. It is a young PNe, 
~1000 years since it left the AGB stage, characterized by a relatively compact ionized region 
(Volk & Kwok 1997) driven by the radiation that arises from a central source at ~200,000 K 
(Latter et al. 2000). In addition, a rich C-chemistry is taking place in the larger molecular 
envelope (~40" in CO; Masson et al. 1985) that surrounds the inner photon dominated 
regions. In fact, NGC 7027 was the first object to show UIB emission latter attributed to 
the PAH emission (Gillctt et al. 1973) and the only source where pure rotational lines of 
CH+ have been detected (Cernicharo et al. 1997). 
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Therefore, the detection of unidentified features in the IR spectrum of NGC 7027 has 
historically contributed to a better understanding of the interstellar carbon complexity. In 
this work, we report the detection of a new UfIB at ~98 /xm (~102 cm~^) in the far-IR 
spectrum of NGC 7027 and speculate about it possible origin. 

2. Observations 

The possible emission/absorption produced by the low-energy bending modes of C„ 
species has been searched using the Long-Wavelength Spectrometer (LWS; Clegg et al. Swin- 
yard ct al. 1996). In the case of the NGC 7027 PNe, we have used all the LWS AOT LOl 
data taken by ISO (sec Herpin et al. 2002). The LWS grating spectra of NGC 7027 taken 
during orbits 21, 342, 349, 356, 363, 377, 537, 552, 559, 566, 579, 587, 594, 601, 706, 713, 
720, 727, 734, 741, 755, 762, 769, 776, and 783 have been averaged. The total on-source 
time was 53409 s and the signal-to-noise (S/N) ratio is very high. Here we present part of 
the spectrum taken with the LWS between ~86 and 110 //m (~116-91 cm~^) at a resolution 
of A/AA ~ 200. The resulting spectrum is shown in Fig. 1, middle panel. We also present 
the same wavelength range in the LWS spectrum [TDT19800158] of IRC+10216 (see Fig. 1, 
lower panel). The full spectrum was shown and modeled by Cernicharo et al. 1996. 

All data were processed following pipeline number 9 and analyzed using the ISO Spectral 
Analysis Package^ (ISAP). Typical routines include: deghtching spikes due to cosmic rays, 
oversampling and averaging individual scans and removing baselines polynomials. 

In addition, we present part of the Short-Wavelength Spectrometer (SWS; de Graauw 
et al. 1996) spectrum of NGC 7027 between ~4-7 //m (~2500-1430 cm"^). The resolution 
of the SWSOl observation template in this interval is A/AA ~ 1500. Many carbon chains 
(polar and non-polar) possess prominent IR active stretching modes in this wavelength range. 
However, Fig. 2 shows that the bulk of the emission in NGC 7027 arises from narrow ionic 
and recombination lines from the inner H II region. Only the well-known UIB broad emission 
at ~6.2 fim is detected. This band corresponds to the relaxation of excited vibrational states 
of aromatic species pumped by absorption of visible-UV photons from the hot central star. 
The full SWS spectrum has been presented by Bernard Salas et al. (2001). 



^The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instruments 
Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC,MPE, RAL, and SRON. 
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3. Results 

After identified and modeled the molecular rotational lines of CO J=30-29 to 25-24, 
CH+ J=4-3 and OH J=5/2-3/2 (see Herpin et al. 2002 for the details), a strong 

band-like feature constituted by several lines was identified between 91 and 102 /xm. Only 
the high-J CO rotational emission produces modest contribution in the band wavelengths. 
Fig. 1 also shows the maximum contribution expected for the rotational 5/2-3/2 fine 

of OH at ~98 //m estimated from the clearer detections at ~79 and ~84 //m. Note that 
no water lines have been detected in the far-IR spectrum of NGC 7027 (Cernicharo et al. 
1997), and thus, no more 0-bearing species are contaminating the ~86-110 /im window. 

Due to the presence of an H ll region in the inner envelope, we have also searched for all 
the possible atomic transitions that could arise from the nebular gas. The most likely atomic 
lines in a PNe such as NGC 7027 (fine structure and recombination lines), taken from the 
Atomic Line List by Peter van Hoof have been analyzed. We note that lines from [Viv], 
[Vlll] and [Tim], with transitions in the considered wavelength range, have been previously 
detected in the optical spectrum of NGC 7027 (Baluteau et al. 1995), while several Hi 
recombination lines are observed in the mid-IR (e.g. see Fig. 2). 

To estimate the Hi emission in the far-IR we have considered a Hll region (2" — 4") 
of ~0.022 MQ(see the model of Volk & Kwok 1997) with nH=ne and assumed that the 
Hi population is given by the Saha-Boltzmann equation. With these parameters, we can 
predict the Hi recombination lines observed in the mid-IR by Bernard Salas et al. (2001), 
and estimate the opacity of the far-IR lines. We found that only the Una series with n 
from 10 to 15 may produce modest emission in the LWS range. However, none of these 
lines appear between 90 and 110 /im. In addition, we have analyzed the possible overlapping 
with fine structure lines. In particular, the [Till]96.68, 106.27 and [Vll]97.79 /im lines could 
produce some emission in the considered range. Even assuming a vanadium abundance 10 
times larger than the solar abundance, the [Vll]97.79 /im (^03-^02) line would be extremely 
weak. Note that the [Vll]141.68 fim (^Da-^Di) line is neither detected in NGC 7027. Finally, 
the solar abundance of Ti is ~10 times that of vanadium. We also computed that is too low 
to produce significant emission in the [Till]96.68, 106.27 lines. 

Hence, the integrated band intensity and the different lines or sub-bands could not be 
fully assigned to any of these atomic lines nor to the pure rotational lines with significant 
line strength arising from the light species in the circumstellar envelope. 

We also note that after subtracting the crowed emission of CO, HCN and HCN^,i;, (pure 
rotational lines in vibrational excited states) from the spectrum of IRC+10216, some of the 
remaining peaks also agreed with those unidentified peaks observed in NGC 7027 (see Fig. 
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1). However, the spectral resolution is too poor to distinguish unidentified features from the 
forest of HCNt,j5 rotational lines that contaminate the far-IR spectrum of IRC+10216 but 
are missing in NGC 7027. The HCN abundance in IRC+10216 is very high, ~3xl0^ with 
HCN/CO=1/10. In fact, HCN is the main coolant of this C-rich AGB star (Cernicharo 
et al. 1999). However, the molecule is photodissociated during the post-AGB evolution as 
the UV radiation field from the evolving star increases. From far-IR observations of HCN 
pure rotational fines, Herpin et al. (2002) found minimum abundances of HCN/CO=1/100 
and 1/1000 for the proto-PNe CRL 2688 and CRL 618 respectively. The extreme case is 
NGC 7027, were a great part of HCN molecules must have been photodissociatd in CN 
and H. In fact, Bachiller et al. (1997) found CN/HCN~10 for this object. Thus, the low 
abundance of HCN compared to IRC+10216 is consistent with the non detection of far-IR 
HCN lines (Herpin et al. 2002). Pure rotational lines of other cyanopolyynes species such as 
HC3N can not produce emission in this range because of the high energy levels associated 
with its far-IR transitions. 

Taking into account the high S/N ratio of the NGC 7027 spectrum, and the absence 
of the band in sources such as Orion, Sgr A or 0-rich evolved stars, we believe that the 
lines are real. Due to the chemistry of NGC 7027, the most probable carrier should be the 
low-energy bending mode of a polyatomic molecule containing carbon. 

4. Discussion 

Both PAHs and pure carbon chains possesses low-energy bending modes in the far-IR. 
Contrary to the UIB emission in the mid-IR, the far-IR skeletal modes of PAHs (vibrations 
associated with the bending of the skeletal structure) depends of the exact nature of the 
species. Hence, the rovibrational structure of the PAH skeletal mode has to be calculated 
for each specific PAH. As a example, the lowest energy active mode of coronene (C24H12) 
fies at ~127 cm~^ (Jobfin et al. 2002), while for ovalene (C32H14), two far-IR skeletal modes 
at ~120 cm~^ and ~66 cm~^ have been theoretically investigated (Mulas et al. 2003). The 
far-IR emission of a few PAHs has also been observed in gas phase experiments (Zhang et 
al. 1996). However, if one of the observed UflBs belongs to a specific PAHs, other bands 
(generally stronger) should be observed in the IR spectrum. This is not the case of the new 
UfIB presented in this work and we have considered the low-energy modes of pure carbon 
chains as possible carriers of the observed UflB. 

The detection of several rovibrational lines of the V2 bending mode of C3 (~63 cm~^) in 
Sgr B2 and IRC+10216 evidenced high abundances for this species and opened the possibility 
to detect more carbon chains in the ISM and GSM (Cernicharo et al. 2000). At the low 
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resolution of the LWS/grating, the more intense C3 hnes [Q(2,4,6)] are blended with the 

strong [C 11] 158 /xm hne and no assignation could be made in NGC 7027. However, the 
tentative detection of the 2/5 bending mode of C4 in many objects including NGC 7027 
(~174 cm~^; Cernicharo et al. 2002), if confirmed, will imply a C3/C4<10 abundance ratio, 
which will suggest that the abundance of C„ does not decrease drastically when n increases. 

In spite of their importance, the C„ bending modes are difficult to characterize spec- 
troscopically in the laboratory or by ab initio computations. As expected, the main uncer- 
tainties difiiculting their astronomical detection are the band origin and the IR intensity. 
Moreover, the induced dipole moments and other spectroscopic constants are often known 
for the fundamental transition but less or nothing is known about overtone transitions within 
the bending mode or transitions to other excited states including the stretching modes. 

The ground electronic state is different for the odd- and even- numbered cumulenic 
chains. This property determines the observed band shape of the bending mode. Odd- 
numbered linear C„'s have a singlet ground electronic state which results in a per- 
pendicular ^Uvib — ^T^vib vibronic spectrum with a strong -Q branch and weak -P and -R 
branches. On the other hand, even-numbered linear C„'s have a triplet ground electronic 
state. Hence, each rotational line is split in three components due to spin couplings. This 
splitting increase with the number of atoms but is particularly small for C4 and Cq (Giesen 
et al. 2001). The vibronic transition is now '''Tlvih — ''''^vih, and despite the null component of 
the electronic orbital momentum in the ground state, the spin-orbit constant (A50) could 
be large in even- C„'s. In such a case, the resulting band-shape will be notoriously affected 
by the value of Aso- 



4.1. Tentative detection of Cg and C5 

With the exception of the z/2 mode of C3, the band origins of other C„ bending modes are 
not accurately constrained. Hence, their assignation as carriers of the UflBs is not obvious. 
Although the low-lying bending modes of C5 and Cq have not been directly observed in gas 
phase, different studies suggest wavenumbers around ~100 cm^^ (see below). The observed 
new UfIB around ~98 //m is composed by several peaks at ~91.8, 92.8, 93.9, 95.2, 97.6, 98.9, 
100.5 and 101.7 /im. In the following, we will consider the possibility that any of these peaks 
is related with the ug and i/^ bending modes of Cq and C5 respectively. 

C5 was first detected in the gas phase by Heath et al. (1989) who observed the 1/3 
stretching mode at ~2169 cm~^ (~4.6 jim). The same mode was observed in IRC-l-10216 by 
Bernath et al. (1989). Actually, C5 is the largest C„ detected in the GSM. On the other hand. 
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Cq was first identified in the laboratory by its electronic spin resonance spectrum in an Ar 
matrix (Van Zee et al. 1987) and later observed in gas phase through its z/4 stretching mode 
at ~1960 cm~^ (~5.1 /xm) by Hwang et al. (1993). Moazzen-Ahmadi et al. (1989) observed 
the (z/3+z/7)-z/7 and (z/3+2z/7)-2z/7 hot bands arising from the bending mode of C5. From 
the /-doubling constant q^, they estimated a frequency of z/7=118±3 cm"-*^. More recently, 
C5 has been observed in photoelectron spectra. In particular, Arnold et al. (1991) measured 
the 21/7 transition and estimated i/7=101±45 cm~^, while Kitsopoulos et al. (1991) obtained 
i/7=107±5 cm~^. Ab initio calculations predict an infrared intensity around 36 km mol"^ 
[A(i/7=l-0)~0.06 s-^; Hutter & Liithi 1994; Martin et al. 1995]. Fig. 1 shows the expected 
band-shape for a ^11^6 — ^"^vib transition with this intensity and the molecular constants for 
the 1^7 mode (Moazzen-Ahmadi et al.) centered at 104.8 cm^^ (95.4 /im). As said before, 
any of the band peaks observed in NGC 7027 around ~98 fim (~102 cm^^) could well be 
responsible of the Uj mode. However, only the best fit to the band including the Cg model 
is shown in Fig. 1. If the carrier of the 95.5 fim feature is finally C5 at z/7=104.8 cm~^, the 
estimated column density assuming an excitation temperature of 100 K (typical of a PNe 
envelope) is iV(C5)=1.8xl0-'^^ cm~^. Assuming 1.5-4.0 magnitudes of visual extinction in 
the neutral envelope (Hasegawa et al. 2000), the typical abundance of C5 in NGC 7027 
would be (0.5-1.0) xlO-^ similar to that derived by Bernath et al. (1989) in IRC+10216. 

Fig. 2 shows that, for the same column densities required to reproduce the bending 
modes, the U4 and stretching modes of Cg and C5 are not detected in NGC 7027. Differ- 
ent excitation conditions of the stretching and bending modes and geometrical effects may 
explain this conjecture. Although carbon clusters have a moderate size, we could also ex- 
pect an efficient pumping of the stretching and bending modes through UV photons. Mid-IR 
photons could also contribute to the excitation of the high-Uhending excited states but also 
of the stretching modes. Finally, the low-energy bending modes could be pumped alone by 
the absorption of far-lR dust photons (see the case of the C3 excitation in Cernicharo et al. 
2000). However, the smaller volume of the photon-dominated regions relative to the far-lR 
dusty envelope in NGC 7027, and the larger mid-IR flux of IRC-M0216 (where the C3 and 
C5 stretching modes were detected) compared to NGC 7027 (a factor ~10^ at 5 //m, but 
only a factor ~4 at 100 //m), can favor the observation of the C„ bending modes within 
the large ISO beam. Under these conditions, the detection of the bending modes and the 
non-detection of the stretching modes in NGC 7027 can be plausible. 

The ug low-energy bending mode of Cg is even less known. Ab inito calculations 

predict a Ug^lOS cm~^ frequency and an infrared intensity of ~25 km mol~"^ [^(^'g = 
1 — 0)~0.035 s~^; Martin et al. 1990; 1995), while photoelectron spectroscopy experiments 
estimate ug=97±A5 cm"^ (Arnold et al. 1991) and z/9=90±50 cm'^ (Xu et al. 1997). Theo- 
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retical predictions have to be taken into account with caution since they rely in the harmonic 
approximation which can be in error for large amplitude benders such as radicals with low- 
lying electronic states. In particular, the infrared intensity has to be considered as a lower 
limit (and column densities as upper limits). Ab initio calculations do predict a low-lying 
electronic state for open-shell species such as C4 or Cg. The term value for the ^A^ lowest 
energy excited state of Cq has been recently determined by Xu et al. only at ~1400 cm~^ 
from the ground. Hence, Aso can effectively be large even in the ground state (A=0). As 
an example, the Aso constant of C4H (^E ground state) is very large, ~3 cm~^ (Yamamoto 
et al. 1987), because its ^11 electronic excited state is only ~468 cm~^ above the ground. 
In the case of even-numbered chains, Aso can easily be larger because of the higher 
spin multiplicity of the ground state and the larger A value of the lowest-lying electronic 
excited state. Fig. 1 also shows the expected band-shape for a ^Umb — ^"^vib transition with 
As'o=3.5 cm^-*^ (we estimated A5o~4 cm~^ for C4), Cg molecular constants from Hwang et 
al. (1993) and van Zee ct al. (1987;1988), and the z/g band origin at 98.3 cm~^ (101.7 fim). 
Assuming the same excitation temperature than for C5 we derive A^(C6)=0.8xlO^^ cm~^. 
Hence, the Cq abundance will be a factor 2 smaller than that of C5. 

Much more work has to be done to fully understand and characterize the low-energy 
vibrations of C„ chains in order to assign the UflBs observed by ISO. The cyclic C„ isomers 
could also produce spectral features in the far-IR, however, their active modes are even 
less known. Waiting for such progresses, space observations offer the unique opportunity 
to obtain spectra in the far-IR domain where the bending modes appear and to motivate 
further laboratory and theoretical studies. The presence of small C„ chains in the space 
as well as their high reactivity, suggest that these species can be involved in the formation 
of more complex organic molecules. In fact, for cumulenic clusters with n=l(} to 20, linear 
structures are thought to close into rings, while for n>30, these species are thought to be 
more stable in aromatic and fuUerene-like structures (O' Brien et al. 1987). 

Future space heterodyne telescopes, such as the Herschel Space Observatory^ with much 
better sensitivity and spectral resolution in the far-IR, should allow the detection of longer 
C„ chains through their low-energy bending modes. This will be the fingerprint needed to 
understand the formation and the nature of the UIBs and UflBs carriers. 
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Fig. 1. — Top panel: Expected band shape for a ^Hvib — ^^vib transition (as the 1/7 mode 
of C5) and for a ^Hyib — ^^vib transition (as the ug mode of Cq) with A5o=3.5 cm~^. The 
excitation temperature for both transitions is 100 K. The intrinsic hne width is 10 km s~^ 
and the spectral resolution A/AA is 200. — Middle panel: Observed ISO/LWS spectrum 
of NGC 7027 between ~86 and 110 fim. The ordinate corresponds to the flux over the 
continuum flux and the abscissa to the wavelength in /xm. Large Velocity Gradient model 
for the pure rotational emission of CO, CH+ and OH (Herpin et al. 2002) and the total 
model for the bending modes of C5 and Ce- The thick horizontal arrow represents the 
observational uncertainty for the C5 band origin (from Kitsopoulos et al. 1991). The 
experimental uncertainty for the ug Cq mode is 90±50 cm~^ (Xu et al. 1997). See text. — 
Bottom panel: Observed ISO/LWS spectrum of IRC+10216 between ■--'86 and 110 fxia. The 
ordinate corresponds to the continuum subtracted flux and the abscissa to the wavelength 
in /im. Large Velocity Gradient model for the pure rotational emission of CO z/=0, 1, ^"^CO 
z/=0, HCN and H^^CN u =0, and HCN^ib z/2=l, 2 and z/i,3=l (Cernicharo et al. 1996). The 
rotational transitions of CO, HCN in the ground state, and HCN^j^ ^2=^, 2 and z^i,3=l are 
indicated. 
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Fig. 2. — Observed ISO/SWS spectrum of NGC 7027 between ~4 and 7 /im. The ordinate 
scale is F/F^ and the wavelength is in nm. The fine structure and recombination lines from 
the ionized region are labeled. The expected S^j^ — T^yHy parallel spectra from the z/3 C3, z/3 
C4, 1^3 C5, z/4 Cq and 1/4 C7 stretching modes are also shown. The excitation temperature 
for all bands is also 100 K. The intrinsic line width is 10 km s~^ and the spectral resolution 
A/AA is 1500. Column densities are A^(C5)=1.8xlO^^ cm~^ (from the z/7 bending mode), 
iV(C6)=0.8xlO^^ cm~2 (from the z/9 bending mode), iV(C3) = 10xiV(C5), iV(C7)=iV(C5)/10 
and A^(C4) =0.8x10^^ cm'^ (from Cernicharo et al. 2002). 



